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ABSTRACT 


Results are presented for an exploratory comparative analysis of radiator geometric 
characteristics based on a simplified radiator model with either a pure gas or a 
graphite-gas suspension working fluid. Radiator variables investigated were gas-flow 
Reynolds number, pres sure -drop fraction, and suspension -particle loading ratio. For 
the pure gases, type of gas and Reynolds number had a pronounced effect on tube and 
panel geometry, but Reynolds number and pressure drop fraction had a relatively small 
effect on planform area for fixed overall cycle loss pressure ratio. Only moderate re- 
ductions in radiator planform area (10 to 17 percent) were indicated with the use of 
graphite -gas suspensions. 
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ANALYSIS OF RADIATOR CHARACTERISTICS OF A SINGLE -LOOP 
100-KILOWATT BRAYTON SPACE POWER SYSTEM USING 
A PURE GAS AND A GAS-SOLID SUSPENSION 
by Robert Pfeffer, * Salvatore Rossetti, * and Seymour Lieblein 

SUMMARY 

An exploratory comparative analysis of radiator characteristics for a 100-kilowatt 
single-loop Brayton space power cycle has been conducted based on a simplified radiator 
model in which either a pure gas or a graphite-gas suspension was used as the working 
fluid. Gases considered were helium, neon, and argon, with parametric variations of 
gas-flow Reynolds number, overall cycle loss pressure ratio, fractional pressure drop 
across the radiator, and suspension-particle loading ratio. Recently developed corre- 
lations were used for the required properties of the gas -solid suspension. Outputs of the 
calculations were optimum cycle temperature ratios, radiator planform area, tube ge- 
ometry (number, length, and inner diameter), panel planform area, and optimum 
suspension-particle loading ratio. In addition, an approximate radiator weight compari- 
son was provided for the suspension cases. 

The principal results of the analysis for the assumptions and conditions covered in- 
dicated that, for the pure-gas working fluids, the gas Reynolds number and radiator 
pressure-drop fraction had a relatively small effect on radiator planform areas as long 
as the overall cycle loss pressure ratio was kept constant. However, even at a constant 
cycle loss pressure ratio, type of gas and Reynolds number had a prounced effect on the 
radiator tube and panel geometry. 

Moderate reductions in radiator planform area compared with the pure-gas case 
were obtained, with the use of a graphite-gas suspension, with an average reduction of 
17 percent for helium, 13. 5 percent for neon, and 10. 5 percent for argon. Somewhat 
larger reductions in weight were indicated by the approximate weight analysis, with 
magnitude depending on the basis of comparison between the pure-gas and suspension 
cases. It was concluded that a more precise evaluation of the potential reduction in ra- 
diator area and weight achievable with the use of a solid suspension would require more 
accurate property data for particle suspensions at low loading ratios and more sophisti- 
cated radiator analyses. 


Associate Professor, City College of the City University of New York. 
^Graduate Student, City College of the City University of New York. 


INTRODUCTION 


Advanced space -power generation systems must be able to continuously supply elec- 
trical power for long periods of time. One such power generating system currently 
under technical evaluation is the indirect conversion, closed- loop heat engine based on 
the Brayton cycle. In this system, heat is generated in a nuclear or solar source and 
rejected by a radiator, with power being obtained from a turbine located in the working 
fluid cycle. 

Several analyses of potential performance of a Brayton cycle using an inert gas as 
the working fluid have been conducted (refs. 1 to 4). These studies have indicated that, 
because of the relatively low temperature levels in the radiator and the inherently low 
heat-transfer coefficients of gases, the Brayton-cycle waste-heat rejection structure in 
a single-loop system is relatively large and heavy. 

One method of decreasing the radiator area, and hence its weight, in a single-loop 
system is to increase the heat-transfer coefficient of the gas working fluid, thus allow- 
ing greater heat transfer per unit of radiation area. For a pure gas working fluid, the 
heat-transfer coefficient is governed by the Reynolds number of the flowing fluid, which, 
in turn, is related to the geometric characteristics of the radiator. 

Decreasing the pressure drop through the heat-transfer components of the cycle is 
still another way to decrease the area and weight of the system. Reductions in pressure 
drop will increase the cycle efficiency, which, in turn, causes a reduction in the re- 
quired radiator area. Here too, the radiator flow Reynolds number is a prime influenc- 
ing factor. 

Another possibility would be to decrease the isentropic specific-heat ratio of the 
working fluid. A decrease in the specific-heat ratio of the fluid results in a higher ther- 
mal efficiency for the Brayton cycle and hence a lower required radiator area (ref. 6). It 
has been shown (ref. 5) that the addition of small particles to a gas will produce a sub- 
stantially lower isentropic heat-capacity ratio than that of the pure-gas component. A 
gas -solid suspension will also increase the heat-transfer coefficient of the flowing fluid 
(ref. 5). Thus, the use of a gas-solid suspension as the working fluid in a single-loop 
Brayton cycle may have a potentially advantageous effect on the cycle heat-rejection 
system. A preliminary analysis of such a cycle, based on suspension-property assump- 
tions and simplified conditions made in reference 6, indicated a large reduction in radi- 
ator area. 

A rigorous parametric study of the effects of such factors as gas -flow Reynolds 
number, pressure loss, and particle addition on the geometric characteristics of a 
Brayton cycle radiator over wide ranges of variation is a cumbersome task in view of 
the complex nature of the design procedures for a gas radiator (e. g. , ref. 7). However, 
if a valid simplified radiator model is used, an analysis procedure can be developed that 
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can provide first-order or preliminary comparative information with only a moderate 
calculation effort. 

This report presents the development and method of solution for a combination cycle 
analysis and radiator model that can be used to determine the preliminary geometric 
characteristics of radiator designs using either a pure gas or a gas-solid suspension as 
the working fluid. The effects of pertinent variables such as gas -flow Reynolds number, 
pressure loss through the radiator, and the properties of the pure gas or gas-solid sus- 
pension on the radiator design are determined and compared over a wide range of vari- 
ables. The calculations are made for the specific case of a representative 100-kilowatt- 
output power cycle with fixed maximum temperature and fixed component efficiencies. 
Argon, neon, and helium are considered for the cycle working fluid, and graphite is con- 
sidered as the suspension particle material. 


THERMODYNAMIC ANALYSIS 

A schematic flow diagram of a single-loop Brayton cycle is shown in figure 1. The 
numbers correspond to the state -point designation used in this analysis. The gas leav- 
ing the heat source at point 1 expands through the turbine to point 2, thereby producing 
the mechanical work necessary to drive the compressor and alternator. From the tur- 
bine the gas enters the recuperator where it is cooled to point 3, as it transfers heat to 



Figure 1. - Flow diagram of single-loop Brayton cycle. (Numbers 
denote state-points.) 
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the gas leaving the compressor. Final cooling to point 4 takes place in the radiator 
where the excess heat is rejected to space. The gas leaving the radiator is then com- 
pressed to point 5, heated in the recuperator to point 6, and further heated back to 
point 1 in the heat source. 

Equations for calculating the thermodynamic cycle efficiency rj , the working fluid 
flow rate W, and the prime radiating area (symbols are defined in appendix A) of 
the Brayton cycle shown in figure 1 have been developed in reference 1 and are pre - 
sented in appendix B. The prime radiating area A^ in equation (B4) is defined as the 
area of either a tubular radiator without fins or a fin-and-tube radiator with 100-percent 
fin efficiency. It is the minimization of this area which is generally used as a criterion 
for optimum cycle temperature selection. The prime area is an acceptable preliminary 
criterion for the evaluation and comparison of radiator characteristics because it can 
be obtained analytically in a simple manner, is proportional to radiator weight and size 
(radiator planform area), and minimizes at approximately the same cycle tempera- 
tures as the full size and weight of a fin-and-tube radiator (ref. 7). 

Examination of the equations of appendix B indicates that, for fixed cycle power 
output P g , maximum temperature Tj, space sink temperature T , and component 
efficiencies t)q, Ty-p, and E, the prime radiating area will be determined by the fluid 
isentropic specific-heat ratio y, the radiator heat-transfer coefficient h^, the cycle 
loss pressure ratio t^/vq, and the temperature ratios T£/Tj and T^/Tj. The 
values of the last two quantities are generally selected according to some cycle de- 
sign criterion (i. e. , minimum prime area), so that they become secondary variables. 
Thus, the primary factors of interest in determining radiator characteristics are 
r T / r c> and y- 

The cycle loss pressure ratio r^,/r^, defined as the ratio of turbine inlet to exit 
pressure divided by the ratio of compressor exit to inlet pressure, is also equal to the 
product of the ratios of exit to inlet pressure for all the heat-transfer components. 
Consequently, represents the fraction of compressor pressure ratio that can be 

recovered to do the work in the turbine and is an indicator of the cumulative pressure 
drops of the heat-transfer components. Prime radiating area is very sensitive to the 
value of the cycle loss pressure ratio (e. g. , ref. 1). Therefore, the effect of the con- 
tribution of the pressure drop through the radiator to the overall loss is an important 
factor for consideration. 

Another important parameter necessary to determine the prime radiating area is 
the gas film or inside-surface convective heat-transfer coefficient h. . This quantity 
is related to the radiator heat-transfer coefficient h R , which is required in computing 
the prime radiating area (eq. (B4)) by the relation 
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( 1 ) 



where A. is the inside surface area of the tube and is the prime radiating area. 
References 1 to 3 assumed h^ to be an independent parameter that was held constant 
in the cycle calculations. It is known, however, that h^ is dependent on gas Reynolds 
number, tube diameter, and fluid properties. This is especially true if a gas-solid 
suspension is used as the working fluid, because h^ will depend strongly on the sus- 
pension properties. For these reasons, h R will not be assumed constant but will be 
computed in the cycle analysis for different conditions and different working fluids. 

The specific heat C p has a unique value for a pure gas at a given temperature. 
However, for a gas-solid suspension, the value of the specific heat will depend on the 
specific heats of the carrier gas and the particle and the particle loading ratio 77 , de- 
fined as the ratio of the mass flow rate of solid particles to that of the pure gas, 


w„ 


The specific heat of the suspension is given by (ref. 5) 



1+677 
I + 77 


( 2 ) 


( 3 ) 


where 6 is defined as the ratio of specific heat of the particles to that of the pure gas 



The specific heat of a gas -solid suspension can either be greater or smaller than the 
value for the carrier gas (ref. 5), but because the product wC^ appears as the variable 
affecting the prime area in equation (B4), the specific heat is not involved as a cycle 
parameter (eq. (B3)). However, the specific heat of the working fluid must be known in 
order to determine the mass flow rate w (eq. (B3)). 

The proper value to be used for the isentropic heat ratio of a suspension v is not 
as clearly defined as for the specific heat. If, however, equilibrium between gas and 
particles can be assumed that is, T = T and V = V , then the isentropic specific- 

o P & 
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heat ratio of a suspension y g in terms of the readily available properties of the system 
can be estimated by (ref. 5) . 



For a suspension, the specific-heat ratio will be less than that of the carrier gas. For 

large suspension-particle loading ratios ( 7 / > 10), y_ approaches 1. 0, and the suspen- 

s 

sion behavior approaches that of an incompressible fluid. 

The effect of a reduced isentropic specific-heat ratio y on the thermodynamics of 
the Brayton cycle is to increase the cycle efficiency and hence reduce the heat rejection 
required by the radiator. The general sensitivity of the prime radiating area Ap to the 
isentropic specific-heat ratio is illustrated in figure 2 for representative cycle condi- 
tions and a fixed value of h R . It should be noted, however, that a reduction in isentropic 
specific-heat ratio results in an increase in the pressure ratios across the turbine and 



Isentropic specific-heat ratio, y 


Figure 2. - Variation of minimum prime radiating area with isentropic specific- 
heat ratio. Inlet temperature, 2200° R (1220 K); turbine efficiency, 0.85; 
compressor efficiency, 0.80; convective heat-transfer coefficient, 10 Btu per 
hour per square foot per °F (56.8 W/(nr)(K); recuperator effectiveness, 

0.85; surface emissivity, 0.86; shaft power, 100 kilowatts. 
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Isentropic specific-heat ratio, y 


Figure 3. - Variation of turbine pressure ratio with isentropic specific-heat ra- 
tio, with turbine temperature ratio as parameter for fixed turbine efficiency 
of 0.85. 


compressor for a given temperature ratio. This effect on the turbomachinery is illus- 
trated in figure 3. 


RADIATOR ANALYSIS 
Radiator Configuration 


The radiator configuration considered in the analysis (fig. 4) is one with a central 
fin-tube panel geometry of axial length L and N number of tubes. This configuration 
is composed of tubes of inside diameter D i and outside diameter D Q and rectangular 
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fins of thickness and width 2 If. Heat rejection occurs from the outside surface of the 
tubes and from both surfaces of the fins by radiation to space. For simplicity, the 
headers are neglected, and only the fin tube panel is considered. 

Preliminary design calculations of the 100-kilowatt single-loop Brayton cycle radi- 
ators based on a comprehensive computer program (ref. 7) had indicated that, for tube 
inside diameters larger than around 3/8 to 1/2 inch (0. 95 to 1. 27 cm), the ratio of fin 
half width to tube outside diameter If/ D q , and the overall radiating effectiveness rj f 
showed little variation with tube inside diameter D. and type of inert gas used for 
minimum -weight configurations. It was also acceptable to consider a constant average 
value of the ratio of tube outside to inside diameter D q /Dj. The use of a constant value 
of Dq/Dj is consistent with internal pressure considerations and, in the more exact 
calculations, showed small effect on the variation of planform area with tube inside di- 
ameter. Thus, it was established that, for comparative analysis purposes, a simplified 
radiator fin tube model can be established based on representative constant values of 
If/ D q , rjf, and D q /D^ that would produce radiator planform areas closely representative 
of those of minimum -weight radiators. 
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Based on this simplified model geometry, equations will now be written that deter- 
mine the radiator pressure relations, heat-transfer rate, and radiator panel geometry. 


Pressure Relations 


In order to be able to investigate the effect of the pressure drop through the radiator 
on the cycle, it is necessary to analyze the pressure drops through all the heat-transfer 
components in relation to the cycle loss pressure ratio. The cycle loss pressure ratio 
and the pressure at various points in the cycle can be expressed in terms of the pressure 
drop across the radiator (AP) rac j. The development of these equations is given in appen- 
dix C. From these derivations, it is shown that, 



( 6 ) 


where Kj is defined as the ratio of the pressure drop in the radiator to the total system 
pressure drop, 


K 


1 “ 


< AP >rad 

< AP >tot 


( 7 ) 


and Kg is defined as the ratio of the pressure drop in the low-pressure side of the re- 
cuperator to the total pressure drop in the cold side of the cycle (radiator and recuper- 
ator), or 


(AP) P 9 - P„ 

K„ s §. (8) 

(AP) c P 2 - P 4 

In the analysis, Kj and Kg are taken as independent parameters, within certain bounds 
(see appendix C), and Pj is considered to be fixed. The turbine pressure ratio Pj/Pg 
is given by 



0) 




-(y/y- 1) 


Equations for the total pressure at all cycle stations in terms of the preceding variables 
are also derived in appendix C. 

The cycle loss pressure ratio can also be expressed in terms of the radiator 
pressure-drop parameters and a fixed value of the total pressure loss in all the cycle 
components except the radiator (AP) com p- This expression (eq. (C18)) is 


(1 - k 2 ) 


a - k., 

K., 


(AP) 


comp | 


(i - k 2 )| 


f 1 - K, 


K, 


(AP) 

j v 7 comp 


- Kr 


Rap) 


( 10 ) 


comp 


The radiator pressure drop can then be related to the radiator geometry, the fluid 
properties, and the cycle mass flow rate (appendix C) as 


(AP) d = 2F[l + 4. 0(Re)"°‘ 32 J|o. 046(Re)" 0 ' 2 V— V— ^ 

g 8 \ D iA g c P 3/ 




( 1 + ^) 4 of N 


(ID 


where 



and 


w 

( Re ) g = 5 

M(1 + V) ~ DjN 
4 1 


( 12 ) 


( 13 ) 
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For the pure gas, rj = 0, and w = w . Thus additional equations are required to relate 

® & 

the geometric parameters of the radiator (N, Dj and L). 


Heat-Transfer Coefficient 


Although previous cycle analyses assumed a constant value for the heat-transfer 
coefficient related to radiating area hjj (eq. (1)), in the calculation of (eq. (B4)), the 
use of the radiator model allows for the calculation of this factor for different gases (or 
gas-solid suspensions), different flow conditions, and different temperatures and pres- 
sures. The heat-transfer coefficient between the gas and the tube inner surface hj, can 
be computed for a pure gas from the relation (ref. 8 ) 


= 0. 023(j5yRe)“- 8 (Pr )°- 3 


w 


where (Re)^ is given by equation (13) for 77 = 0 and w g = w^, and 

( c p) "g 
(Pr) = SIlsJ. 


(14) 


(15) 


Reference 5 has shown that the relation between h , the heat-transfer coefficient 

s 

of a gas-solid suspension, and h , the heat-transfer coefficient of a pure gas, can be 

o 

taken as 




-0.32 

g 



(16) 


so that, for the case of a suspension flow, 



( 17 ) 


The values of h., using either equation (13) for the pure-gas equation, or equa- 
tion (17) for the gas -solid suspension, are converted to h R by equation ( 1 ) where the 



■ ■ ■■ 
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tube inside surface area is 


Aj = jtDjLN 


(18) 


The prime radiating area A is obtained directly from the cycle analysis as indicated 

XT 

in equation (B4). It now remains to relate the calculated prime area to the geometry of 
the radiator model. 

The total radiating surface area Ap of the radiator model (fig. 4) is given by 


A 


R 



(19) 


The prime area, which is the effective radiating area, will differ from the total radiating 
area as given by equation (19) because of two factors: (1) The fins are not radiating at 
100-percent efficiency, and (2) some of the radiation leaving the radiator surface is in- 
tercepted by some other part of the surface (i. e. , the view factor to space is less than 
one). 

The first factor can be accounted for by correcting the fin radiating area by a factor 
less than 1. 0 (the fin efficiency). Reference 9 has shown that the second factor can be 
corrected by using the projected area of the tubes 2D q LN rather than the actual radiating 
area 3 rD Q NL. The prime (or effective) radiating area can therefore be described by 


Ap = 2D 0 NL ♦ 



( 20 ) 


so that, in terms of D^, 


A 


P 



Therefore, from equations (18) and (21) 



( 21 ) 


( 22 ) 
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so that equation (1) becomes, 



D iV * D o/ 


( 23 ) 


The use of either equation (14) or (17) for h. requires an expression for the tube 
inside diameter Dj. This relation is obtained by combining equations (11) and (13) to 
give 


°i “ l/- Re - E - 0. 092(Re)~°- 2 [l + 4. 0(Re)~°' 32 1 ,] — ” ^ 
f ap rad 8 8 J S c p 3 


An expression for the length of tubes L required in equation (24) can be obtained by 
combining equations (13) and (21) to yield, 


L = 


A p (Re >g 

8w g 


Mff( 1 + 7?) 



(25) 


Thus, if representative values of rjp ^/D q , and D Q /D i can be prescribed, equa- 
tion (23) allows hj^ to be related to gas properties and the radiator geometry. 


Panel Geometry 

The use of the simplified radiator model allows the calculation of significant features 
of the panel geometry such as the planform area, panel aspect ratio, and area vulnerable 
to meteoroid damage as functions of the tube flow Reynolds number, the fluid properties, 
and the cycle conditions. 

The planform area is defined as the projected area of the radiator panel neglecting 
the headers, and is described by 


A , = NLD 
pi o 



(26) 
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where the number of tubes N is then obtained from equation (C40) as, 


N = 


4w_ 


ffG g Df(l + V ) 


The planform area is related to the prime area by the ratio, 


( 27 ) 


[ 1+2 — 

V_ 1 D - 


Ap 2 f 1 + 2t >r^] 


(28) 


The vulnerable area (area vulnerable to damaging impact by meteoroids) is 


A, = ttD NL 
v o 


(29) 


so that in ratio form 


pl 


1 + 2 — — I 


D 


(30) 


Another significant factor considered in the design of radiators is the panel aspect 
ratio <p which is defined as the length-to-width ratio of the panel and is expressed as, 


(p = 


N (D 0 + 2 ^) A pl 


(31) 


CALCULATIONS 


Method of Solution 


Once the cycle parameters r? T , E, and i* T /r c are specified and the independent 
variables Tg/T^, T^/T^, P^, and P g are chosen, the cycle is fully described for a 
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given pure gas or a specified gas-solid suspension. From the cycle analysis, w, A p , 
and all temperatures are determined. When the radiator model constants (D 0 /D., 
rjj, Kj, and Kg) are fixed and the gas Reynolds number in the tubes of the radiator are 
prescribed, the four equations (13), (24), (25) and (27) in conjunction with the cycle out- 
puts were used to determine the four parameters, D^, G^, L, and N, respectively. An 
iterative procedure was required based on an initial assumption of due to the depend- 
ence of hp (which is necessary for the cycle analysis) on D^. All calculations were 
performed on an IBM 7094 computer. The computational procedure is presented in 
appendix D. 


Inputs 

The example chosen for the parametric analysis was a 100-kilowatt shaft power out- 
put cycle using either pure helium, neon, or argon, or a suspension of graphite in these 
gases as the working fluid. From turbomachinery considerations, compatible values of 
2200° R (1222 K) and 25 200 pounds per square foot (12. 1x10^ N/m^) chosen, respec- 
tively, as the turbine entering temperature T^ and pressure P^. The equivalent sink 
temperature was chosen as 400° R (222 K), and the value of 0. 86 was chosen for the 
emissivity, e, of the tubes and fins. Values of the compressor, turbine, and recuper- 
ator efficiencies were taken as 0. 85 ( 77 q = rjrj, = E = 0. 85) for both the pure gases and the 
gas -solid suspensions. 

The three important cycle parameters are Tg/T^, T^/Tj, and the cycle loss pres- 
sure ratio, r T /r c . In this analysis, Tg/T^ was varied between 0. 70 and 0. 90, and 
T^/Tj was varied between 0. 25 and 0. 45 in increments of 0. 025. The allowable or de- 
sign value of the cycle loss pressure ratio is generally a compromise between cycle ef- 
ficiency and component flow passage size and weight. The best or optimum value can be 
determined from a complete system weight and performance analysis; it is not clear, 
beforehand, what an acceptable value might be. Much will depend on which components 
are most sensitive to size. In most cases, the principal components that are dependent 
on flow passage size are the heat source (e. g. , reactor and shielding, radioisotope and 
shielding) and the radiator. In order to ascertain the effect of r^/r^ on the radiator 
design, the cycle loss pressure ratio was taken as 0. 85, 0. 90, and 0. 95 for the pure-gas 
analysis; whereas an average value of 0. 90 was chosen for the gas-solid suspension 
studies. 

As a parametric variable, the Reynolds number was chosen between 10 000 and 
500 000 depending on the gas. The Reynolds number range for each gas was chosen to 
keep the tube inside diameter above 3/8 inch (0. 95 cm) in order not to invalidate the 
model assumptions discussed previously. Based on the analyses of gas-cycle radiators, 


representative geometric parameter values were selected as D q /D^ = 1. 33, l^/. D q = 3. 5, 
and = 0. 67. 

A constant value of 0. 25 was considered appropriate for the recuperator pressure- 
drop fraction Kg. In order to test the sensitivity of the radiator design to the radiator 
pressure-drop fraction Kj, three values of K^ (0. 2, 0. 5, and 0. 75) were selected as 
input data for the pure gas studies. The value of K^ = 0. 75 is the limiting value of this 
parameter when Kg = 0. 25 if no pressure losses occur in either the high-pressure side 
of the recuperator or the heat source. For the gas-solid suspensions studies, a single 
value of Kj = 0. 5 was chosen. 

The significance of the factor as a direct measure of the pressure drop through 
the radiator is illustrated in figures 5 and 6. Figure 5 shows the relation between the 
pressure drop through the radiator ((AP) rac j/Pg) an d Kj for three different values of 
r T / r c f° r the case of Kg = 0. 25. The calculation was made for a value of Tg/Tj = 0. 8. 
It is seen that if the overall cycle loss pressure ratio is fixed, decreasing the value of 
Kj will decrease the pressure drop in the radiator, but will also simultaneously increase 
the pressure drop in the other heat-transfer components. 

In certain applications, it may be necessary to fix the pressure drop in the other 
heat -transfer components rather than the overall cycle loss pressure ratio. Figure 6 



Figure 5. - Variation of ratio of pressure drop 
across radiator to inlet pressure with radi- 
ator pressure-drop fraction for fixed cycle 
loss pressure ratio. Recuperator pressure- 
drop fraction, 0.25; outlet-inlet turbine 
temperature ratio, 0.8. 



Figure 6. - Variation of cycle loss pressure ratio with 
radiator pressure-drop fraction. Recuperator pressure- 
drop fraction, 0.25; outlet-inlet turbine temperature 
ratio, 0.8. 
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shows the relation between the cycle pressure ratio and for Kg = 0. 25 for five dif- 
ferent values of the total pressure drop in the other heat-transfer components, 
(AP) com p/Pj (eq. (10)). Thus, if the total pressure drop in the other heat-transfer com- 
ponents is specified (with Kg = 0. 25), figure 6 can be used to choose a value of K^ 
which will result in an acceptable value of the overall cycle loss pressure ratio. If, on 
the other hand, all three individual pressure drops comprising ( A P) cornp are specified, 
the values of and Kg can be calculated directly from equations (7) and (8), and the 
operating cycle loss pressure ratio is obtained from equation (10). 

The pure gases that were considered are the inert gases: helium, neon, and argon. 
The properties of these gases required as inputs to the analysis were thermal conduc- 
tivity k, viscosity p, specific heat C , isentropic specific heat ratio y, and gas con- 
stant R* . Values of these properties for the temperature range of interest are presented 
in appendix E. The gas -solid suspensions studied were suspensions of graphite in these 
same inert gases. The pertinent physical properties of graphite are also presented in 
appendix E. 


RESULTS AND DISCUSSION 

The geometric parameters of the radiator were computed for each gas and suspen- 
sion case as a function of the temperature ratios Tg/T^ and T^/Tj. Radiator plan- 
form area was used as the principal output factor, and the minimum value of the plan- 
form area was selected as the criterion for the optimum cycle temperature ratio. Once 
the optimum temperature ratios were established, a parametric analysis of the impor- 
tant variables affecting the radiator geometry at these temperatures was made possible. 
These variables include the inert gas to be used as the working fluid, the suspension- 
particle loading ratio, the gas Reynolds number, the cycle loss pressure ratio, and the 
radiator pressure -drop fraction Kj. 

Although radiator planform area is taken as a principal measure of radiator per- 
formance in view of its direct relation to radiator weight and size, it should be noted 
that variations in prime and vulnerable area can also be determined from applications of 
equations (28) and (31). 


Pure-Gas Results 

Opt imum temperature ratios . - A typical plot of the planform area Apj against 
Tg/T^ with T 4 /Ti as a parameter is presented in figure 7 for neon. Although a mini- 
mum value of temperature ratio is indicated in all cases, the plot shows that the curves 
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Temperature ratio across turbine, T 2 /T 1 

Figure 7. - Variation of radiator planform area with temperature ratio across 
turbine for 100-kilowatt neon cycle. Gas Reynolds number, 50 000; cycle 
loss pressure ratio, 0.90; radiator pressure-drop fraction, 0.5; recuperator 
pressure-drop fraction, 0.25. 

are fairly shallow in the vicinity of the minimum. In general, the optimum temperature 
conditions (Tg/T^ and T^/T j) were obtained in the analysis by scanning the computer 
output for the smallest value of V and not by plotting the planform area. The incre- 
ment used for varying Tg/T^ and T^/Tj in the computer program was quite small 
(±0. 025) so that the maximum error from this procedure would be ±0. 025 in the tem- 
perature ratios for the cases where the curves were flat in the vicinity of the minimum. 

A comparison of the optimum temperature ratios Tg/T^ and T^/T^ obtained by 
this work and by the method of reference 1 based on minimum prime radiating area and 
fixed heat-transfer coefficient h R is plotted in figures 8 and 9, respectively, as a func- 
tion of cycle loss pressure ratio. Both calculations used the same temperature ratio 
increments and selection process. The specific values of optimum T2/TJ and T^/Tj 
for helium, neon, and argon fell within the maximum error band, so that a single curve 
can be plotted for all three gases. As can be seen from these figures, optimization of 
temperatures using the radiator model yields slightly higher values of both Tg/T^ and 
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Optimum temperature ratio across 
turbine, T 2 /Tx 



Cycle loss pressure ratio, rj/r^ 

Figure 8. - Variation of optimum temperature ratio across 
turbine with cycle loss pressure ratio. Average for 
helium, argon, and neon. Maximum error band, +0.025 
on temperature ratio. 



Figure 9. - Variation of optimum temperature ratios of 
compressor to turbine inlet with cycle loss pressure 
ratio. Average for helium, argon, and neon. Maxi- 
mum error band, +0.025 on temperature ratio. 



Cycle loss pressure ratio, r T /r c 

Figure 10. - Variation of prime radiating area with cycle loss 
pressure ratio at optimum temperature ratios. Gas Reynolds 
number, 50 000. 
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T^/Tj for all three gases. This is attributed to the difference in the value of h R used 
in each method. A further indication of the effect of the value of h R as computed by 
the two approaches is shown by figure 10. Figure 10 is a plot of prime radiating area 
as a function of cycle loss pressure ratio with Kj as a parameter for argon and helium, 
and of prime radiating area obtained by using a constant value of h R “ 5 and 10 Btu per 
hour per square foot per °F (28. 4 and 56. 8 W/(m^) (K)) included as a comparison. The 
lower values of prime radiating area obtained in this study (for helium and argon) as 
shown in figure 10 are due to the higher values of the actual heat-transfer coefficient 
h R obtained in the analysis compared with the constant values. 

The numerical results also indicated that, within the optimization selection process 
used, neither radiator pressure-drop fraction K^, nor gas Reynolds number has any 
effect on the optimum temperature ratios for a given gas and cycle loss pressure ratio. 
However, as shown in figures 9 and 10, the optimum temperature ratios tended to in- 
crease with increasing cycle loss pressure ratio for all three gases. 



Figure 11. - Variation of radiator planform area with gas Reynolds number. Radiator pressure-drop fraction, 0.5; 
recuperator pressure-drop fraction, 0.25. 
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In order to present smooth curves for the radiator geometry factors (L, N, Dj, 
etc.), some of the operating conditions chosen were not exactly at the precise minimum 
value of the planform area. However, these points were always within ±0. 025 of the 
temperature ratios at the precise minimum and always within 0. 5 percent of the 

value of minimum planform area. 

Planform area . - The effect of gas Reynolds number on the radiator planform area 
at the optimum temperature ratios and fixed radiator pressure-drop fraction K 1 is 
shown in figure 11. As can be seen from this figure, the effect of gas Reynolds number 
is to increase the value of A j monotonically as (Re) g is increased over the range of 
interest for each gas and cyc?e loss pressure ratio considered. The magnitute of this 
increase, however, is small (only around 5 to 10 percent over the range of Reynolds 
numbers considered). 

The variation of radiator planform area with radiator pressure-drop fraction, Kj, 
at a fixed Reynolds number of 50 000 is shown in figure 12 for several values of cycle 



Radiator pressure-drop fraction, Kj 


(a) Helium. (b) Neon. (c) Argon. 

Figure 12. - Variation of radiator planform area with radiator pressure-drop fraction at optimum temperature ratios. Gas 
Reynolds number, 50 000; recuperator pressure-drop fraction, 0.25. 
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loss pressure ratio and component-loss factor ((AP) CQ /P^) for helium, neon and 
argon respectively. The figure indicates that an increase in Kj from 0. 2 to 0. 75 at 
constant cycle loss pressure ratio does not affect the planform area appreciably (less 
than a 10 percent decrease in for the range of covered). However, if a fixed 
component loss is prescribed, a substantial reduction in planform area results from re- 
ductions in because of the effect on r T /r^,. Thus, the planform area is much more 
sensitive to the value of overall cycle loss pressure ratio than to the value of K^. 

Tube geometry . - The variation of tube length with gas Reynolds number at optimum 
temperature ratio conditions and fixed Kj is shown in figure 13 for three gases and 
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Figure 13. - Variation of tube length with gas Reynolds number. Radiator pressure- 
drop fraction, 0.5; recuperator pressure-drop fraction, 0.25. 


cycle loss pressure ratios. The figure shows that L varies linearly with (Re)^, as in- 
dicated by equation (25) . It appears that the use of helium results in much longer tubes 
than the use of argon or neon for the same conditions. The effect of the cycle loss pres- 
sure ratio appears to depend on the gas used; the effect being largest with helium. In all 
cases the tube length decreases with increasing cycle loss pressure ratio. 

The effect of radiator pressure -drop fraction K 1 on the tube length was relatively 
minor (less than 10-percent decrease in length for an increase of Kj from 0. 2 to 0. 75 
compared with the effect of the other factors involved over the range of variables cov- 
ered. 
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Tube inside diameter, 





.2 .4 .6 .8 

Radiator pressure-drop fraction, Kj 


Figure 15. - Variation of tube inside 
diameter with radiator pressure-drop 
fraction. Gas Reynolds number, 

50 000; recuperator pressure-drop 
fraction, 0.25. 
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The gas Reynolds number also has a very strong effect on the tube inside diameter, 
as illustrated in figure 14. The linear dependency of D. on Reynolds number can be seen 
from equation (24) because L varies linearly with Reynolds number. The figure also indi- 
cates the large effect that the gas exerts on the tube diameter; the tube diameters needed 
with helium as the working fluid are almost four times those required using argon as the 
working fluid for the same Reynolds number and cycle loss pressure ratio. 

The effect of the radiator pressure-drop fraction Kj on the tube inside diameter is 
shown in figure 15, with r^/r^, as a parameter at a constant (Re)^ of 50 000. As Kj 
decreases (corresponding to decreasing values of the pressure drop across the radiator) 
the tube diameters become larger as would be expected for a lower pressure drop. The 
effect of the cycle loss pressure ratio, appears generally small in comparison with the 
effects of Reynolds number and gas used. 

Equation (27) shows that the number of tubes is inversely proportional to the square 
of the tube diameter. Because the tube diameter is almost linear with Reynolds number, 
N will vary approximately with the inverse square of the Reynolds number. Because 
either decreasing r^/r^, or increasing results in a higher pressure drop across 
the radiator (see fig. 3), the number of tubes will tend to increase with a decrease in 



Figure 16. - Variation of number of radiator tubes with gas Reynolds number. 
Recuperator pressure-drop fraction, 0.25. 
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r T / r c or an increase in K^. However, the variation of N with t^/yq and was 
relatively small. Figure 16 illustrates the effect of a combination of these two variables 
that results in the maximum variation in the number of tubes. It is thus seen that the 
two major factors affecting the tube geometry (L, and N) are the gas used and the 
gas flow Reynolds number. 

Heat-transfer coefficient. - The major factors which influence the heat-transfer 
coefficient based on radiating area h^ appear to be the gas used, the radiator 
pressure -drop fraction Kj, and the gas Reynolds number. Plots of h R against (Re)g 
for helium, neon, and argon, are shown in figure 17 for a fixed value of r^/r^. Al- 
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(b) Neon. 


Figure 17. - Variation of heat-transfer coefficient with gas Reynolds number. 
Cycle loss pressure ratio, 0.90; recuperator pressure -drop fraction, 0.25. 
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though hp is proportional to (Re)^' , it is also inversely proportional to D^, which 
increases "with increasing (Re)g ( see fig- 14). The net effect is that h R decreases with 
(Re) g as indicated in the figure. The fact that higher values of resulted in higher 
heat-transfer coefficients was expected because a high value of K ^ represents a large 
allowable pressure drop through the radiator tubes and hence results in a small tube 
diameter. Similar trends of variation of h^ were obtained at the other values of the 
cycle loss pressure ratio. 
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Figure 18. - Variation of panel aspect ratio with gas 
Reynolds number. Radiator pressure-drop fraction, 0. 5. 
Recuperator pressure-drop fraction, 0.25. 


Panel aspect ratio. - Figure 18 shows the variation of the panel aspect ratio 
(eq. (31)) against gas Reynolds number with the cycle loss pressure ratio at a fixed 
value of K r If radiator design considerations limit the aspect ratio over a certain 
range (e. g. , 0. 5 to 2) then these curves indicate the appropriate Reynolds number range 
of operation. The figures show that the major influences on the aspect ratio are the gas 
used and the Reynolds number. At a given Reynolds number, the use of helium results 
in a panel aspect ratio approximately 10 times that of neon and 100 times that of argon. 
Increasing the cycle loss pressure ratio (at constant (Re)^) decreases the aspect ratio 
as indicated in the figures, whereas varying Kj had a negligible effect on the aspect 
ratio (less than 10 percent for a change in Kj from 0. 2 to 0. 75). Thus, even at low 
radiator pressure drops (high r T /r c and low K,), reasonable panel aspect ratios 
(0.5 to 2) are obtainable for all three bases if the Reynolds number is appropriately 
chosen. 

Gas velocity. - The effect of gas Reynolds number and radiator pressqre-drop frac- 
tion on the average and inlet velocity through the tubes is illustrated in figures 19 
and 20, respectively. The inlet velocity is approximately 20 percent higher than the 
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Figure 19. - Variation of velocity of gas in tubes with gas Reynolds number. Cycle loss pressure 
ratio, 0.90; radiator pressure-drop fraction, 0.5; recuperator pressure-drop fraction, 0.25. 



Figure 20. - Variation of average velocity 
with radiator pressure-drop fraction. 
Gas Reynolds number, 50 000; re- 
cuperator pressure-drop fraction, 0. 25. 



average velocity. The velocity along the tubes using helium is more than double that 
using argon, all other factors being equal. Figure 19 shows that the velocity increases 
slightly with Reynolds number. Figure 20, shows that a decrease in from 0. 75 
to 0. 2 results in a decrease in the average velocity because the allowable pressure drop 
in the tubes is reduced; whereas, an increase in r^/r^, from 0. 85 to 0. 95 causes a 
decrease in velocity for similar reasons. Because all the velocities computed were 
substantially below the acoustic velocity, the assumption of incompressible flow used 
in the analysis is acceptable. 


Gas-Solid Suspension Results 

The effect on radiator characteristics of the use of a gas-solid suspension as the 
cycle working fluid was determined for three suspensions: graphite -helium, graphite- 
neon, and graphite-argon. Inasmuch as the primary emphasis in these calculations was 
on the effect of the particle loading ratio of the suspension and since previous results 
indicated that the effects of r T /r c and were generally smaller than those of the 
gas used and the gas -flow Reynolds number, the calculations were conducted for a single 
value of cycle loss pressure ratio of 0. 90 and radiator pressure-drop fraction of 0. 5. 

The minimum value of the planform area was selected as the criterion for determining 
both the optimum cycle temperature ratios and the optimum particle loading ratio for each 
gas. As indicated previously, property data for the suspensions were obtained from the 
correlations of reference 5. 

Optimum temperature ra tios. - Plots of optimum Tg/T^ and T^/T^ against par- 
ticle loading ratios are shown in figures 21 and 22. As can be seen in the figures, both 
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Figure 22. - Variation of optimum ratio of compressor inlet to turbine inlet temperature with 
particle loading ratio. Cycle loss pressure ratio, 0.90; radiator pressure-drop fraction, 0.5; 
recuperator pressure-drop fraction, 0.25. 


of the optimum temperature ratios increase with increasing particle concentration. The 
rate of increase is smallest with helium and most pronounced with argon. The increase 
in the optimum temperature ratios with increasing particle loading ratio is due to the 
fact that the isentropic specific-heat ratio y_ decreases with particle loading ratio. 
Therefore, the temperature difference across both the turbine and compressor decrease. 
The disparity in the optimum temperature ratios for the three different suspension 
carrier gases is due to the difference in the specific heat of the gases resulting in dif- 
ferent values of 6 for each suspension. 

Planform area . - The variation of panel planform area with particle loading ratio at 
several gas Reynolds numbers for three gas-graphite suspensions is plotted in figure 23. 
Values of the planform area for 77 = 0 represent the pure -gas case. 

The figures show that a minimum planform area is obtained at a specific particle 
loading ratio for each gas, and this optimum particle loading ratio varies only slightly 
with the gas Reynolds number. For the helium suspension, the optimum particle loading 
ratio is between 4. 5 and 6 . 0 for the range of Reynolds numbers covered. For the neon 
suspension, the optimum particle loading ratio is around 0. 75 and 1. 0, and, for the 
argon suspension, the optimum particle loading ratio falls between 0. 30 and 0. 40 for the 
respective ranges of representative gas Reynolds number. A plot of optimum particle 
loading ratio against Reynolds number is shown in figure 24 for all three gases. 

A minimum of the planform area with increasing particle loading ratio occurs be- 
cause of the opposing effects of the isentropic specific-heat ratio and the convective 
heat-transfer coefficient h^. For a given gas, an increase in the particle loading ratio 

decreases y , increases the cycle efficiency, and thus tends to decrease the planform 
s 

area. At the same time, an increase in 77 decreases h R (as will be subsequently ex- 
plained) and, therefore, tends to increase the planform area. 

Because y and h„ are both functions only of the product of the particle loading 
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(c) Argon -graphite suspension. 

Figure 23. - Variation of planform area with particle loading ratio for different gas Reynolds numbers. Cycle loss pressure 
ratio, 0.90; radiator pressure-drop fraction, 0.5; recuperator pressure-drop fraction, 0.25. 
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Figure 24. - Variation of optimum particle loading ratio with gas Reynolds 
number. Cycle loss pressure ratio, 0.90; radiator pressure-drop fraction, 
0.5; recuperator pressure-drop fraction, 0. 25. 
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Figure 25. - Variation of planform area with 
product 5 t]. Gas Reynolds number, 

50 000; cycle loss pressure ratio, 0.90; 
radiator pressure-drop fraction, 0.5; 
recuperator pressure-drop fraction, 0.25. 


co 



ratio times the ratio of the heat capacities of the solid to the gas 677, it is the value of 
this quantity which determines the minimum planform area. This is illustrated in fig- 
ure 25, which shows that, at a gas Reynolds number of 50 000 for all three graphite-gas 
suspensions considered, the minimum planform area occurs at approximately the same 
value of 677 equal to 1.4. 

The magnitude of the maximum decrease in panel planform area potentially achiev- 
able with the use of a gas -solid suspension is more graphically illustrated in the plots of 
the ratio of planform area of the suspension to that of the pure gas A ^(77 = opt)/A p j(7j = 0) 
(fig. 26). Compared with the pure- gas, average reductions in panel planform area from 
17 percent for helium to 10 percent for argon are indicated for the range of variables covered. 



Figure 26. - Variation of ratio of planform area at optimum particle loading ratio to planform area at 
zero particle loading ratio with gas Reynolds number. Cycle loss pressure ratio, 0.90; radiator 
pressure-drop fraction, 0.5; recuperator pressure-drop fraction, 0.25. 


The empirical friction factor and heat-transfer correlations used in the suspension 
radiator calculations contained extrapolated values for loading ratios lower than 1. 0 as 
discussed in reference 5. Therefore, there is some uncertainty in the absolute magni- 
tude of the resultant values of radiator geometry and in the optimum values of particle 
loading ratios as determined herein. 

Tube geometry. - The effect of particle loading ratio on tube geometry is illustrated 
in figures 27 to 29 by plots of the calculated variations of N, D^, and L at representa- 
tive optimum particle loading ratio against gas Reynolds number as compared with the 
corresponding values for the pure-gas case (77 = 0). The results in all cases are for the 
optimum values of ^ 2^1 and ^ 4^1 for each suspension system. It. is seen from 
figure 27, that the use of an optimum particle loading ratio markedly reduces the re- 
quired number of tubes at a given Reynolds number. The ratio between the suspension 
and pure-gas cases averages between 0. 25 and 0. 36 for the range of conditions covered. 

Figure 28 shows an increase in tube length with the use of an optimum, loading ratio. 
The average ratio of suspension length to pure-gas length varied from around 1. 55 to 
1. 75. The suspension also results in an increase in tube inside diameter for a given gas 
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Figure 27. - Variation of number of radiator tubes with 
gas Reynolds numbers at representative particle loading 
ratios. Cycle loss pressure ratio, 0.90; radiator 
pressure-drop fraction, 0.5; recuperator pressure-drop 
fraction, 0.25. 
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Figure 28. - Variation of radiator tube length with gas 
Reynolds number at representative particle loading 
ratios. Cycle loss pressure ratio, 0.90; radiator 
pressure-drop fraction, 0.5; recuperator pressure- 
drop fraction, 0. 25. 





Figure 29. - Variation of radiator tube inside diameter 
with gas Reynolds number at representative particle 
loading ratios. Cycle loss pressure ratio, 0.90; radi- 
ator pressure-drop fraction, 0.5; recuperator pressure- 
drop fraction, 0. 25. 
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Reynolds number (fig. 29) . The ratios in this case also varied somewhat with the gas 
used: around 1. 95 for helium, 1. 60 for argon, and 1. 75 for neon. 

Panel aspect ratio . - The panel aspect ratio is also greater at a given Reynolds 
number with the use of a gas -solid suspension (see fig. 30). Average ratios vary from 
around 2. 7 for argon to 3. 5 for neon. The increase in panel aspect ratio is primarily 
the result of the reduction in number of tubes and increase in tube length as indicated in 
figures 28 and 29. If the radiator panel aspect ratio is restricted to a certain range of 
values because of the specific application of the radiator, figure 12 can be used to deter- 
mine the approximate Reynolds number range of operation. For example, if the panel 
aspect ratio is required to be 2. 0, the operating gas Reynolds number should be around 
11 000 for helium-graphite, 33 000 for neon-graphite, and 105 000 for argon-graphite 
suspensions at optimum conditions. 

Heat-transfer coefficient. - The effect of the gas -solid suspension on the heat- 


transfer coefficient based on radiator area h R is illustrated in figure 31. The values of 
h_ for the suspension are lower than for the pure gas. This apparent anomoly is due to 
the relative variations of Ir for the pure gas and for the suspension. Although, the 


ratio h g /h increases for a given tube diameter as loading ratio is increased (as indi- 
cated by eq. (16)), for the optimum loading ratios considered, the tube diameter is 


greater than for the pure-gas case at a given Reynolds number (fig. 29). Thus, there is 


a net decrease in the value of h„ compared with 

b 


h 


g 


at a given gas flow Reynolds 
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Figure 30. - Variation of panel aspect ratio with Reynolds 
number for representative particle loading ratios. 

Cycle loss pressure ratio, 0.90; radiator pressure-drop 
fraction, 0.3; recuperator pressure-drop fraction, 0.23. 
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(a) Helium and argon. 
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Gas Reynolds number, {Re) g 


(b) Neon. 

Figure 31. - Variation of radiator heat -transfer coefficient with Reynolds number at 
representative particle loading ratios. Cycle loss pressure ratio, 0.90; radiator 
pressure-drop fraction, 0.5; recuperator pressure-drop fraction, 0.25. 


number, and, with Aj/A^ a constant in all cases (eq. (22)), there is a net reduction in 
hj^ for the suspension compared with the pure-gas case. 

Weight comparison . - The results of figure 26 indicated a relatively modest re- 
duction in radiator planform area with the use of a graphite-gas suspension at a repre- 
sentative optimum particle loading ratio. Although the radiator model and computation 
procedure were not originally devised to investigate radiator weight, it is possible, on 
the basis of the present analysis, to obtain a comparison of approximate weights of pure- 
gas and suspension radiators. The development of the weight equations for the radiator 
fin tube model of figure 4 is detailed in appendix F. From this development, it is shown 
that the ratio of suspension radiator weight is 



in which values of A^ and are obtained from the comparable calculated results of 
the pure-gas and suspension cases (i.e. , figs. 11, 23, 26, and 29). Thus, for the pre- 
scribed values of (p o /D.) and Z f /D Q , and the assumption that P f /p t , Z f /D Q are the 

same for the pure-gas and suspension cases, ratios of radiator weight can be obtained 

from equation (32) if the value of (d/d. ) /(l) /D.) is known or prescribed. 

s g 

Although Dq/D. has been specified to be 1. 33, for both the pure-gas and suspen- 
sion, it is likely that this value might be reduced for the case of the suspension because 
a larger tube inside diameter is obtained at any fixed Reynolds number (fig. 29). Thus, 
it would be of interest to plot the variation of weight ratio W g /W g as a function of 

the ratio (d q /d^ ^(D q /D. ) for several representative values of fin thickness parameter 

tf/D . Calculations of these variations for the three gases were made for 

= 1.33, Z f /D 0 = 3.5, Pf/p t - 1-0 and 3.4, and t f /D Q = 0.02 and 0.04. Com- 
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(b) Neon; optimum suspension -particle loading ratio, 1. 



(c) Argon; optimum suspension-particle loading ratio, 0.375. 


Figure 32. - Variation of radiator weight ratio with outside-to-inside tube diameter ratio. Gas 
outside-to-inside tube diameter ratio, 1.33; fin length to tube outside diameter ratio, 3.5; 
fin length to tube outside diameter ratio, fin thickness to tube outside diameter ratio, and 
fin to tube density ratio are same for gas and suspension. 
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parisons were then made on the basis of the same gas Reynolds number for the pure -gas 
and suspension cases, and also on the basis of the same panel aspect ratio for the pure- 
gas and suspension cases, as shown in figure 32. Also shown in the figure are the 
values of tube diameter ratio for the criterion of equal meteoroid protection for the 
pure-gas and suspension cases. These values represent likely or realistic values of 
suspension tube diameter ratio (d /D.) as obtained from consideration of meteoroid 

protection theory in determining tube wall thickness (appendix F). 

Figure 32 shows that a substantial reduction in radiator weight is possible with the 
use of a suspension working fluid if the tube wall thickness can be decreased to give re- 
duced values of (d q /d. ) . For the practical case of tube wall thickness required for 

s 

equal meteoroid protection, it appears that the weight reduction for the suspension radi- 
ators for constant aspect ratios is somewhat greater than for the constant Reynolds 
number. For the fixed Reynolds number cases, the weight reduction varied from around 
2 to 20 percent (while the planform area reduction (fig. 26) varied from around 10 to 
18 percent). For a panel aspect ratio of 2. 0, the weight reduction for the suspension 
radiators compared with the pure-gas radiators was around 15 to 25 percent (compara- 
ble planform area reduction is 12 to 19 percent). For both comparison bases, helium 
showed the largest weight reduction percentage, and argon showed the smallest reduc- 
tion. 

In view of the reduced values of D q /D. indicated for the fixed Reynolds -number 
suspension radiators when compared for equal meteoroid protection (around 1 . 16 to 
1. 20 compared with 1. 33 for the pure-gas radiator model), check calculations were made 
to determine the effect of a reduced value of D q /D^ on the calculated planform area of 
the suspension radiators. These results indicated an average reduction in suspension 
radiator planform area of around 2 percent. A more precise determination of planform 
area and weight reductions for the suspension radiators will require a more exact radia- 
tor analysis such as that illustrated by reference 7. 

Compari son witji results of reference 6 . - An analysis of the use of a suspension of 
graphite particles in argon, neon, and helium as the working fluid in a Brayton space 
power cycle was also conducted in reference 6 . The analysis of reference 6 determined 
prime radiating area. For the cycle and assumptions considered, the analysis showed 
first that the prime radiating area decreased rapidly as the product 6 77 was increased 
to about 3. 0 and then decreased much more slowly with 677 . The results obtained in this 
report indicated that, at a Reynolds number of 50 000, the planform area showed a mini- 
mun at a 677 of 1. 4. The second result from reference 6 was an indicated reduction in 
prime radiating area (at 677 = 3. 0) of around 67 percent compared with the pure-gas 
cases. Inasmuch as the current analysis revealed only a 10 to 18 percent reduction in 
planform area at optimum 677 , it is important to explore the basis for the differences. 
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The analysis of reference 6 was based on the equations for cycle efficiency and 
prime radiating area developed in reference 1 and presented in appendix B. Equilibrium 
conditions were assumed to exist in all cycle components, and the isentropic specific- 
heat ratio of the suspension was the same as equation (5). However, the suspension 
heat-transfer coefficient h R used in the equation for prime radiating area was taken to 
vary linearly with 677, according to 


“’’Mg <33) 

For the analysis of this report, in contrast, it was found that the radiator heat-transfer 
coefficient h R , based on the suspension data correlations of reference 5, actually de- 
crease with particle loading ratio (fig. 31) because of the increase in tube inside diam- 
eter (eq. (17)). Thus, lower values of prime radiating area showing no minimum with 
loading ratio can be expected from the analysis of reference 6. 

A second difference between the two cycle analyses is the use, in reference 6, of a 
recuperator scheme in which two zero-gravity gas-solid separators were used. One 
was located at the outlet of the turbine to bypass the hot particles directly to the inlet of 
the heat source. The second separator was located at the outlet of the compressor to 
bypass the cooled particles back to the radiator. Thus, the particle heat recovery oc- 
curred at 100-percent efficiency. Inasmuch as the particles constitute a large fraction 
of the mass of the working fluid, the overall effectiveness of the recuperator is in- 
creased. In using this scheme, the overall recuperator effectiveness with the suspen- 
sion becomes 


e s = =±ia 04) 

S 1+677 

where E is the normal recuperator effectiveness. Thus, the overall recuperator effec- 
tiveness increases asymptotically toward one as the product of 677 is increased. 

In order to determine the effect on radiating area of the increased recuperator ef- 
fectiveness resulting from the use of the particle separator scheme, a set of calcula- 
tions were performed with the combined radiator and cycle analysis developed in this 

report, and with the recuperator effectiveness E defined by equation (34) . Results for 

s 

a helium- graphite suspension at a particle loading ratio of 7 and a gas Reynolds number 
of 10 000 showed a decrease in planform area of around 32 percent over the pure-gas 
case compared with the 18-percent decrease obtained without particle separators. 

Thus, a large portion of the approximately 67-percent reduction in prime radiating area 
quoted in reference 6 is probably due to the higher values of hp produced by the use of 
equation (33). It appears, therefore, that an accurate evaluation of the relative effec- 
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tiveness of the use of a gas -solid suspension working fluid will require an accurate radi- 
ator analysis method as well as accurate property data for particle suspensions in the 
loading-ratio region of interest. 


SUMMARY OF RESULTS 

An exploratory comparative analysis of radiator characteristics for a 100 -kilowatt 
single -loop Brayton space power cycle has been conducted based on a simplified radiator 
model in which either a pure gas or a graphite -gas suspension was used as the working 
fluid. Results of the calculations for pure helium, neon, and argon at optimum cycle 
temperature ratios showed principally that: 

1. The radiator planform area was determined primarily by the magnitude of the 
overall cycle loss pressure ratio and the type of gas used. The fractional pressure drop 
across the radiator had a significant effect on radiator planform area only if the overall 
cycle loss pressure ratio was varied as well. For a given value of overall cycle loss 
pressure ratio, the radiator flow Reynolds number and fractional pressure drop had rel- 
atively small influence on the radiator planform area. 

2. The radiator gas -flow Reynolds number had a pronounced effect on the radiator 
geometric characteristics of panel aspect ratio and tube inside diamter, number, and 
length. For a given Reynolds number, the values of the geometric characteristics 
varied substantially among the three gases considered. For a fixed panel aspect ratio, 
the required flow Reynolds number was greatest for argon and least for helium. 

For the graphite -gas suspension radiators, principal results indicated that: 

3. An optimum particle loading ratio corresponding to minimum planform area was 
found which varied with the gas used and the flow Reynolds number. Values of optimum 
cycle temperature ratios at these conditions were greater than for the pure-gas case. 

4. Moderate reductions in planform area compared with the pure-gas case were ob- 
tained, with an average reduction of 17 percent for helium, 13. 5 percent for neon, and 
10. 5 percent for argon over the range of conditions covered in the analysis. 

5. An approximate weight evaluation showed somewhat larger percentage reductions 
with the use of a graphite suspension, depending on the basis of comparison employed 
(e. g. , fixed Reynolds number, fixed panel aspect ratio, etc.). For a design panel aspect 
ratio of 2. 0, weight reductions of around 22 percent for helium, 20 percent for neon, and 
15 percent for argon were indicated for the condition of equal meteoroid protection. 

Inasmuch as the results of this exploratory analysis depend on the accuracy of the 
correlations for the pertinent properties of the gas-solid suspension and also of the ge- 
ometric properties of the radiator model, a more precise comparative evaluation will 
require more detailed and sophisticated radiator and cycle analysis methods as well as 
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more accurate information concerning the properties of particle suspensions at low par- 
ticle loading ratios in the optimum region. Finally, even if the more precise analysis 
continues to reveal a significant potential gain in the radiator area and weight with the 
use of particle suspensions, the ultimate feasibility of such an approach requires con- 
sideration of the practical problems involved in the use of such systems. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 15, 1968, 

120-27-04-63-22. 
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APPENDIX A 


SYMBOLS 

2 2 

A- radiator tube inside surface area, ft ; m 

2 2 

A prime radiating area, ft ; m 

^ 2 2 
A , radiator planform area, ft ; m 

2 2 

Aj^ total radiator outside surface area, ft ; m 

2 2 

A y radiator vulnerable area, ft ; m 

C fluid heat capacity, Btu/(lb mass)(°R); W/(m 2 )(K) 

.H 

Dj radiator tube inside diameter, ft; m 

D radiator tube outside diameter, ft; m 

E recuperator effectiveness; defined by eq. (B2) 

F correction factor for nonisothermal flow in radiator 

f Fanning friction factor 

f^ friction factor for a gas -solid suspension based on gas density; defined by 

eq. (C25) 

2 2 

G mass velocity, lb mass/(sec)(ft ); kg/(sec)(m ) 

2 2 

g c gravitational constant, (lb mass/lb force) (ft/sec ); (kg)(m)/(N)(sec ) 

h convective heat-transfer coefficient, Btu/ (ft 2 ) (°R) ( hr) ; (W/(m 2 )(K) 

h| heat-transfer coefficient based on tube inside area, Btu/(ft <s )(°R)(hr); 

W/(m 2 )(K) 

9 n 

hp heat-transfer coefficient based on radiator prime area, Btu/(ft )(°R)(hr); 

WAm 2 )® 

Kj radiator pressure-drop fraction: ratio of pressure drop through radiator to 

total pressure drop in cycle 

K2 recuperator pressure-drop fraction: ratio of pressure drop through recuper- 

ator to cycle cold side pressure drop 

k thermal conductivity of gas, Btu/(sec)(°R)(ft); W/(m)(K) 

L length of radiator tubes, ft; m 

half width of fin, ft; m 
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N 

P 

Pr 

P s 

(AP) C 

^^comp 

(AP) h 


(AP) 

(AP) 

(AP) 


rad 

rec 

tot 


Q rad 

R' 

Re 

r T / r c 

T 


w 


L t 

V 

W 

w 


y 

6 

e 

V 


number of radiator tubes 

O 

gas pressure, lb/sq ft; N/m 
Prandtl number 
output shaft power, kW 

sum of pressure drops across radiator and low pressure side of recuper- 

p 

ator (cold-side pressure drop), lb/sq ft; N/m 

sum of pressure drops of heat-transfer components, excluding radiator, 
lb/sq ft; N/m 2 

sum of pressure drops across heat source and high pressure side of recu- 

p 

perator (hot-side pressure drop), lb/sq ft; N/m 

2 

pressure drop across radiator, lb/sq ft; N/m 

2 

pressure drop across cold side of recuperator, lb/sq ft; N/m 

total pressure drop through heat-transfer components in the cycle, lb/sq ft; 
N/m 2 

total heat rejected from radiator 

gas constant, (lb force) (ft)/(lb mass) (°R); j/(kg)(K) 

Reynolds number 

cycle loss pressure ratio, P 1/^2/^ 5/^4 

absolute temperature, °R; K 

equivalent sink temperature, °R; K 

absolute temperature of tube outer wall, °R; K 

thickness of fin, ft; m 

thickness of radiator tube, ft; m 

velocity through radiator tubes, ft/sec; m/sec 

radiator weight, lb mass; kg 

mass flow rate, lb mass/sec; kg/sec 

isentropic specific-heat ratio 

ratio of particle to gas specific heats 

emissivity of radiator tubes and fins 

loading ratio of suspension, mass flow rate of solids/mass flow rate of gas 
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compressor efficiency 
thermodynamic cycle efficiency 
fin efficiency 
turbine efficiency 

gas viscosity, lb mass/ (ft) (sec); N/(sec)(m 2 ) 

o o 

density, lb mass/ft ; kg/m 

3 3 

density of fin metal, lb mass/ft ; kg/m 

o o 

density of tube metal, lb mass/ft ; kg/m 
Stefan- Bo ltzman constant, Btu/(ft 2 )(°R^); J/(m 2 )(K^) 
radiator panel aspect ratio, ratio of overall length to width 
Subscripts: 

av average conditions 

g gas 

p particles 

s suspension 

1 turbine inlet 

2 turbine outlet 

3 radiator inlet 

4 radiator outlet 

5 compressor outlet 

6 heat-source inlet 


% 

v cy 

V f 

p 

pf 

Pt 

a 

<P 
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APPENDIX B 


CYCLE THERMODYNAMIC RELATIONS 

The thermodynamic relations associated with the Brayton cycle are developed in 
reference 1. The developed relations which are pertinent to the present work are pre- 
sented as follows: 



where the recuperator effectiveness E is given by 



(Bl) 


(B2) 


(B3) 
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where the wall temperature T w at any point along the tube is related to the gas tem- 
perature T by 


T = T w + r 1 ( T w' T *) (B5) 

n R 

A numerical technique is required in the determination of the local wall temperature. 

The total heat rejected from the radiator is equal to the heat lost by the gas flowing 
through it and is given by 

«*rad = » C p (T 3 ‘ T 4> < B6) 
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APPENDIX C 


PRESSURE RELATIONS 


Fixed Cycle Loss Pressure Ratio 

The cycle loss pressure ratio i‘ T /r ( -, is defined as (ref. 1) 



However, 

P 5 = P 1 + < AP >h 


(Cl) 


(C2) 


and 


p 4 = p 2 " ( A P) C (C3) 

where (AP) h is the hot-side pressure loss or the pressure loss through the recuperator 
and heat source, and (AP) c is the cold-side pressure loss or the pressure drop through 
the radiator and recuperator. Therefore, 



(C4) 


Because the turbine inlet pressure P ^ will be assumed fixed in this analysis, equa- 
tion (C4) can be written as, 
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(C5) 


r C 



To enable the pressure drop through the radiator to be related to the cycle loss 
pressure ratio, the pressure drop in the radiator will be expressed as some fraction of 
the total pressure drop through the heat-transfer components in the cycle by the equation 


„ . < AP >rad 

1 < AP W 


(C6) 


In addition, the pressure through the low-pressure side of the recuperator, (AI?) rec will 
be expressed as some constant fraction of the total cold-side pressure drop (AP) c , such 
that 


(AP) rec ( AP ) r ec 

K 9 = — = _ (C7) 

< AP >e <AP) rec + < AP > ra d 

The parameters Kj and Kg as defined are not totally independent of each another. 
Once the value of Kg is fixed, K^ has an upper bound. This upper bound is obtained 
from the assumption of no pressure loss through either the high-pressure side of the re- 
cuperator or the heat source, and is found to be 

K l = 1 - Kg (C8) 

If it is assumed that the pressure losses through both the high- and low-pressure sides 
of the recuperator are the same, then the upper bound on K^ is further reduced and is 
given by 


2K 9 

K n = 1 — (C9) 

1 i + k 2 


based on no pressure losses in the heat source. 
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Since, by definition 


(AP) C 


= < aP >rec + 


(AP) 


rad 


(CIO) 


equation (C7) becomes 


(AP) 


rec 


K 2 

1 - Kr 


(Cll) 


Using equations (CIO) and (Cll) and the equation 

(AP) tot = (AP) h + (AP) C 

yields a rewritten form of equation (C6) 

(AP), 


'rad 


v = K, 


(AP) h + (AP) ra J 


Kr 


1 - Kr 


+ 1 


(C12) 


(C13) 


or 




K, 


1 - Kr 


(C14) 


It now remains to relate (AP) rad to the cycle loss pressure ratio of equation (C5). 
Substituting equations (C7), (Cll), and (C14) into equation (C5) results in 


ll - K r 


(AP) 


rad /^1 


(l _ K 2 A 

^ A ^rad 

V K 1 !-K 2 j 

. P 1 


(C15) 


Rearranging and simplifying equation (C15) yield 
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1 


|(1 - K 2 )P 1 


(C16) 


(AP) 


rad 


r T 

U - K 2 - K x )" 

r C 

K i j 


where Pj is fixed, P 1 /P 2 , the turbine pressure, is given (ref. 1) as 

1 -iy/r- 1 ) 


P l_ 

1 / T 2\" 
l - — (i - — ] 

P 2 

_ V T \ T i/_ 


(C17) 


r x/ r c is prescribed, and Kj and K 2 can be parametrically varied within the limita- 
tion imposed by their definitions. 


Fixed Component Pressure Loss 

The radiator pressure-drop fraction, can also be related to the cycle loss pres- 
sure ratio r T /r£ in terms of a fixed value of pressure loss for all the other heat- 
transfer components, (AP) comp . Because (AP) comp is equal to [(AP) tot - (AP) rad ] by 
using equation (C6) 


^^comp 



rad 


(C17) 


If a fixed value of K 2 is also prescribed, the combination of equations (C13) and 
(C17) gives 


rrp 

r C 



(AP) 


comp 


K, 


(1 - k 2 ) P 2 1 - Kj 


K r 


1 - K„ 


- 1 


2 

(AP) 


comp 1^ 


1 - K. 


and in terms of (AP) comp /P 1 
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r T 

r C 


(i - k 2 ) 


rl - K, 


(AP) 


comp 


Ki 


(i - k 2 ) 


A - K, 
K 1 


l + (AP) comp 


- Kr 


Rap) 


comp 


(C 18) 


Thus, if (AP) com p is specified, the contribution of the radiator pressure-drop variation 
Kj to the cycle loss pressure ratio, and therefore to the cycle efficiency, can be evalu- 
ated from equation (C18). The maximum value of the cycle loss pressure ratio for a 
fixed value of (AP) com p/Pi can be determined from equation (C18) by allowing to 
approach 0 and can be shown to be 


r T _ 1 

r C ! + (AP) comp 


(C19) 


Local Pressures 


If the turbine inlet pressure P^ is prescribed, the pressure at every other point in 
the cycle can be expressed in terms of (AP) rac p Kj, and Kg. 

The pressure at turbine outlet Pg is given by equation (C15). The pressure at the 
radiator inlet is given by 


P 3 = P 2 



rad 


(C20) 


At the inlet to the compressor, the pressure is 


P 4 = P 3 - (AP) ra(J 


and, at the recuperator high-pressure inlet, the pressure is 


1 K 2 

P 5 = Pl + |X. 2 


*1 !- K 2 


- 1 V AP ) rad 


(C21) 


(C22) 
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Finally, if it is assumed that the pressure drop through high-pressure side of the recu- 
perator is the same as through the low-pressure side, the pressure at the inlet of the 
heat source can be expressed as 


P 6= P 1 + 


(1 - K 2 ) - Kj - K x K 2 

K“(l - k 2 ) 


(AP) 


rad 


(C23) 


For the limiting case of no pressure losses through either the high-pressure side of the 
recuperator or through the heat source, Pg = Pg = Pj. 


Pressure Drop Equation 

The pressure drop in the radiator consists of losses through the headers and the ra- 
diator tubes. When it is assumed that the pressure drop through the headers is small in 
comparison with that through the radiator tubes and if entrance, bend, and fitting losses 
in the radiator are neglected, the pressure drop can be expressed by the Fanning equa- 
tion. Thus, 


< iP >rad-»g < CM > 

Because the temperatures and pressure vary along the radiator tube, and f 
must be taken as average values to obtain (AP) rad . The product P g V g however, is 
a constant along the tube. For a gas-solid suspension, equation (C24) is modified to give 

< iP >rad = 2P s(^k 2 ^ < c25 > 


where f* is the effective friction factor of the suspension and V_ is velocity of the 
s s 

suspension. Since in dilute suspensions the volume of solids is small compared with the 

volume of the gas, it can be assumed that V„ can be taken to be the same as the velocity 

s 

of the gas phase V . Furthermore, the suspension-friction factor can be expressed as 

D 



(C 26 ) 
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where f is evaluated at the same gas -state conditions (V, P, T, Dj) as the suspension. 
Thus, equation (C25) becomes, 


< ap >rad 



(C27) 


The friction-factor ratio f' /f„ 

s' g 


is assumed to be adequately represented by the relation 


5 

f. 


-S = 1 + 4. 0(Re) 


-0.32 


g 


V 


(C28) 


presented in reference 5. 

For fully developed turbulent flow in smooth pipes (e. g. , ref. 8), 

f g = 0. 046(Re)’°- 2 


where 


m g 


D p V D.G 
i H g g _ i g 


P 


g 




(C29) 


(C30) 


If it can be assumed that the gas flow is essentially incompressible (there is a ques- 
tion as to what the effective sonic velocity of a suspension is), then the gas phase static 
density is related to the local gas total pressure and temperature by the equation 


P 


g 


P 

R’T 


(C31) 


and the velocity at any point in the radiator tube is given by 


V 


g 


G R’T 

g 

P 


(C32) 


Because equation (C27) requires an effective average velocity, equation (C32) can be ex- 
pressed in terms of the inlet velocity (P = Pg and T = Tg) multiplied by a correction 
factor F derived in reference 7 as 
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so that the effective average velocity is 



FG R’To 
g 3 

P 3 


(C34) 


Thus equation (C27) becomes 


( AP )rad =2F [ 1 + 4 - °( Re )g°' 32 ,][0- 046(Re)~°‘ 2 ] -t 

®i 


G ^'T 3 

«c P 3 


(C35) 


The cycle calculations determine the total mass flow rate w which, in terms of the 
radiator model, is given for a pure gas by 


w g = '’g v g; D ? N (C36 > 

and for a suspension by 

w s - Ps v s l D ? N < C37 > 


For low values of the particle loading ratio (77 < 10) the density of the suspension p is 

s 

approximately equal to (see ref. 5) 

P s = P g (l + 7) (C38) 

With the assumption that V = V , equation (C37) can be written as 

s g 



so that G in equation (C35) can be obtained as 
& 


G 


g 


w s 

(1 + 77)^D 2 N 
4 


(C40) 


Combining equations (C40) and (C35) yields 


(AP) rad =2F [ 1 + 4 - °( Re )g°‘ 32 vjo. 046(Re)'°- 2 ]/^j 


/ R,T 3\ 

w s 

\ S c P 3/ 

(1 + V) 1 D 2 N 


12 


(C41) 


where (Re) is given by, 

o 


w 

(Re) g = (C42) 

M(1 + 7?) -D.N 
4 1 

The radiator pressure drop can thus be related to the radiator geometry, the radiator 
flow Reynolds number, and the cycle mass flow rate. 
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APPENDIX D 


COMPUTATIONAL PROCEDURE 

All cycle and radiator model computations were performed on the IBM 7094 com- 
puter at the Lewis Research Center. The procedure used for the computation is outlined 
herein, 

(1) The gas Reynolds number is set (this number can be varied parametrically). 

(2) The cycle loss pressure ratio r-p/r^ is set (this number can also be varied 
parametrically if desired) . 

(3) The properties of the working fluid are set by feeding in the appropriate data 
from appendix E (depending upon whether a pure gas or a gas -graphite suspension is 
used). If a suspension is considered, the values of r] and 6 are varied parametrically. 

(4) The pressure drop across the radiator is calculated using equation (C14) for the 
specified value of the cycle loss pressure ratio and the prescribed values of Kj, Kg, 
and Pj, and P^/Pg obtained from equation (C15). The pressures at points 3 to 6 in the 
cycle (fig. 1) are calculated using equations (C20) to (C23), respectively. For the limit- 
ing case of Kj given by equation (C8), Pg and Pg are taken to be the same as P^. 

(5) A value of the inside diameter of the tubes is assumed, and h^ is computed 
from equation (23). 

(6) Knowing h^, r^/r^, and the specified cycle inputs as outlined previously, the 
cycle efficiency r] c y, the mass rate of flow w, and the radiator prime area Ap are 
computed (for each Tg/Tj and T^/Tj considered) using equations (Bl), (B3), and 
(B4), respectively. The temperatures at every point in the cycle are also determined. 

(7) Because (Re)^ is specified and Ap and w have been computed in step (6), 
equation (25) allows for the direct evaluation of L (the tube length), once l^/ D Q , D Q /Dp 
and Tfy are specified. 

(8) The value of L is known from equation (25); therefore, D i can be calculated 
using equation (24). 

(9) The value of computed in step (8) is compared with the value of assumed 
in step (5). If the two values differ by more than 1 percent, steps (5) to (8) are re- 
peated using the new value of D. as the assumed value of in step (5). 

(10) When the new value of Dj is close enough to the previous trial to satisfy the 
specifications of the program, the gas mass velocity G is computed from equation (C30) 

o 

and the number of tubes N is computed using equation (26). 

After convergence of D i has been achieved, the computer program prints out the 
prime area Ap, the inside diameter of the tubes Dp the length of the tubes L, the 
number of tubes N, the pressure drop across the radiator (AP) rad , the mass flux of the 
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gas in each G , the tube Reynolds number of the gas (Re) , the heat-transfer coefficient 

B B 

hp, the planform area App the vulnerable area A v , and the average velocity through 
the tubes V ay . The program also prints out all cycle variables and input parameters 
including Tp Tg/Tp T g , 4’ ^^’p^'^s’ 

Tg/Tp T w 3 , Tg/Tp Tg/Tp ?7 cy , and A p /P g , as well as the pressures at each point 
in the cycle’ Pp Pg, Pg, Pp Pg, and Pg, and the pressure ratio across the turbine 

V P 2- 
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APPENDIX E 


PROPERTIES OF GASES 

The values and equations for the gas properties used as inputs to the computer pro- 
gram are listed herein. The equations were obtained from faired curves based on the 
data of reference 10 in the temperature range from 760° to 1960° R (422 to 1090 K). 

The viscosity and thermal conductivity values were determined at the arithmetic mean 
temperature in the radiator. 


Property 


C 

y 


p 


R’ 


k 


c 

y 


p 


r t 


k 


c 

y 


p 


R' 

k 


U. S. customary units 

Helium 


1. 25 Btu/(lb mass)(°R) 

1.67 

386. 25 (lb force/lb mass)(ft/°R) 
0.433 T°- 650 X10' 5 Btu/(sec)(ft)(°R) 
0. 0238 T°- 650 X10' 5 lb mass/ (ft) (sec) 


Neon 


0. 25 Btu/(lb mass)(°R) 

1.67 

77. 25 (lb force/lb mass)(ft/°R) 

0. 0128 T°- 662 xl0' 5 Btu/(sec)(ft)(°R) 
0. 0232 T°- 665 xl0' 5 lb mass/(ft)(sec) 


Argon 

0. 125 Btu/ (lb mass)(°R) 

1.67 

38. 63 (lb force/lb mass)(ft/°R) 
0.00298 T 0-725 xl0' 5 Btu/(sec)(ft)(°R) 
0. 0156 T 0, 665 xl0' 5 lb mass/(ft)(sec) 


SI units 


5230 J/(kg)(K) 

1.67 

2080 j/(kg)(K) 

0. 33 T 0, 650 xl0' 2 j/(sec)(m)(K) 
0. 0519 T°- 650 X10" 5 (N)(sec)/m 2 


1046 J/(kg)(K) 

1.67 

416 J/(kg)(K) 

0. 0098 T 0, 662 X10‘ 2 j/(sec)(m)(K) 
0. 0509 T°- 665 X10" 5 (N)(sec)/m 2 


523 J/(kg)(K) 

1.67 

208 J/(kg)(K) 

0. 237 T°‘ 725 X10" 5 J/(sec)(m)(lO 
0.0355 T 0, 725 X10' 5 CN)(sec)/m 2 
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The only property of graphite used as an input to the computer program is its speci- 
fic heat. This value over the temperature range 500° to 2000° R (278 to 1112 K) was 
taken as 0. 39 Btu per pound mass per °R (1635 J/(kg)(K)). Therefore, the values of 6, 
which were used, are: 


Suspension 

Particle to gas 


specific heat ratio, 


6 

He lium - graphite 

0.312 

Neon -graphite 

1.56 

Argon- graphite 

3. 12 


It was also assumed that the average diameter of the graphite particles was of the order 
of 30 micrometers. 


APPENDIX F 


RADIATOR WEIGHT 


For the radiator fin-tube model of figure 5, the weight of the radiator panels com- 
posed of fins and tubes is given by 



or in terms of ratios 



Now, the planform area is given by 



(FI) 


(F2) 


(F3) 


Substituting in equation (F2) yields 


C i— — i ""x 



Thus, with p t and Z f /D Q taken as the same for both the suspension and pure-gas radi- 
ators, the weight ratio becomes 


w 

w 


S 

g 




(F5) 
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I 


In equation (F5), it is assumed that the parameters p f /p t , f f / D q , and t £ /D Q remain 
the same for the suspension and pure-gas cases. Furthermore, for the prescribed radi- 
ator model, (p o / DjJ = 1.33 and Z f / D Q = 3.5. Thus, with values of (a j {. A and 

/ \ / ✓ \ S S P g 

(Djj /\Pj) obtained from the previous calculation results (e. g. , figs. 24, 25, 27, 
s 'g 

and 30), values of the weight ratio can be obtained as a function of the ratio (D 0 /Dj) f 
( D o/ D il for parametric values of t £ /D Q and p £ /p^ for any desired basis of compari- 
son (e. g. , same Reynolds number, same panel aspect ratio, etc.). Two values of p f /p t 
equal to 0. 34 and 1. 0 were chosen to represent either fins and tubes of the same material 
or aluminum fins with stainless-steel tubes. 

Realistic values of (D /d. ) can be obtained from the calculation results for the 

criterion of equal meteoroid protection. According to meteoroid protection theory, the 
required tube wall thickness is given by 


x. Constant 

(y t ) 1/2 


V 


v0. 25 


-In P(Oh 


(F6) 


where E £ is the modulus of elasticity of the tube material, p £ is the density of the tube 
material, A y is the vulnerable area of the radiator panel, r is the design lifetime of the 
radiator, and P(0) £ is the probability of no puncture during the lifetime of the radiator. 
The tube wall thickness can also be expressed as 



(F7) 


so that taking the ratio of suspension-case tube thickness to pure-gas-case tube thickness 
(from eqs. (F6) and (F7)) results in 



where E £ , p £ , r, and P(0) t are the same for both cases. 
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The vulnerable area of the radiator is given by 

A y = jtD 0 LN (F9) 

From the definition of radiator planform area (eq. (F3)), 



(F10) 


Substitution of equation (F10) into equation (F8) with l f /D Q taken the same in both cases, 
yields the expression for tube diameter ratio for the suspension case for equal meteoroid 
protection 



(Fll) 


In equation (Fll), (d.) /(p ) and (a /( A nl) are obtained from the calculated results 
for the particular basis of comparison desired, while ^D 0 /Dj j = 1. 33 for the pure-gas 


model. 
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